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ABSTRACT: 

In  this  paper  are  presented  the  results  of  numerical  computations  of 
the  effect  of  space  charge  upon  the  image  size  and  image  distance  in  a 
converging  beam  of  charged  particles.    The  minimum  beam  radius  (image 
size)  and  its  position  (image  distance)  are  plotted  as  a  function  of  beam 
perveance  ranging  from  10         to  10      .    The  beam  convergence  angle  was 
varied  from  about  0.03  to  0.  15  radians  and  calculations  were  carried  out 
for  3  values  of  aperture  radius  (0.  1  cm,  0.  25  cm  and  0.  50  cm) .    The 
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I.     INTRODUCTION 

In  recent  years  there  has  been  an  increasing  demand  for  low  energy 
ion  beams,  well  defined  in  space  and  energy,  to  be  used  in  experimental 
investigations  of  atomic  collision  processes.    In  the  energy  region  below 
100  eV  the  effects  of  space  charge  upon  usable  beam  intensity,  beam  shape 
and  ion  energy  become  so  serious  as  to  represent  the  limiting  factors  in 
many  experiments . 

Pierce  has  designed  a  widely  used  lens  utilizing  curved  electrodes 
to  create  the  fields  required  to  simulate  flow  between  concentric  spheres . 
Although  the  Pierce  gun  has  been  used  for  the  production  of  low  energy 
ion  beams  with  some  success,  it  was  recently  shown  by  Simpson  and 
Kuyatt     that  there  exists  a  fundamental  limitation  on  the  beam  intensity 
which  can  be  obtained  by  any  unipotential  gun  operating  at  low  energies. 
They  show  that  in  order  to  attain  the  current  densities  required  for 
saturation,  in  the  presence  of  thermal  effects  at  the  emitter  surface, 
one  would  have  to  exceed  the  space  charge  limitations  imposed  by  the 
diode  nature  of  any  unipotential  gun.    From  this  it  is  concluded  that  the 
range  of  applicability  of  unipotential  guns  lies  at  comparatively  high 
energies  and  wide  beam  convergence  angles.    To  operate  outside  these 
limits  it  is  necessary  to  decouple  the  thermal  energy  and  space  charge 
effects  by  first  accelerating  the  particles  to  a  fairly  high  energy  and 
then  decelerating  them  to  their  final  low  energy.     Simpson  and  Kuyatt^ 
have  designed  an  electron  gun  using  this  multistaging  principle  and  we 
have  adapted  their  design  for  use  with  low  energy  ion  beams. 

A  schematic  diagram  of  the  gun  is  given  in  Fig.  1 .    The  first  stage 
extracts  the  particles  at  high  energies  and  forms  a  beam  which  is  then 
decelerated  to  the  desired  final  energy  in  the  second  stage.    The  first, 
or  extraction  stage  is  based  on  designs  given  by  Soa    and  the  decel- 
eration stage  is  based  on  the  properties  of  plane  electrodes,  as  given 
by  Spangenberg  and  Field. °    The  electrode  potentials  V    and  Va  are 
measured  with  respect  to  the  emitter  which  is  maintained  at  the  desired 
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FIGURE     1 


beam  potential,  V     ,  while  the  exit  electrodes  are  grounded, 
o 

In  addition  to  circumventing  the  basic  problem  described  above,  this 
design  has  considerable  mechanical  advantages  over  the  conventional 
Pierce  gun.    All  electrodes  are  plane;  displacement  of  the  deceleration 
stage  with  respect  to  the  extraction  stage  and  variation  of  the  electrode 
potentials  allows  one  to  obtain  a  wide  variety  of  beam  shapes,  intensities 
and  energies  quickly  and  conveniently. 

The  details  and  results  of  our  experimental  investigation  of  this  gun 
are  reported  elsewhere     and  will  only  be  summarized  very  briefly  here.    To 

obtain  a  low  energy  ion  beam  we  used  a  thermal  emitter  (Li-O'Al   0    *n  SiCL) 

+  7  —8  ~6  9 

as  a  source  of  Li     ions    .    Current  densities  of  10      to  10      amps/cm    at 

beam  energies  of  2  to  100  eV  were  obtained  and  both  the  spatial  and  energy 

distributions  of  the  ions  were  investigated.     Minimum  beam  radii  of  1  to  2 

mm  occurred  a  few  centimeters  beyond  the  last  element  of  the  gun.     Beam 

profiles  were  measured  as  a  function  of  axial  distance  and  well  collimated 

beams  of  up  to  16cm  were  obtained.    The  energy  spread  in  the  beam  was 

found  to  be  about  0.22  eV  in  good  agreement  with  the  expected  thermal 


energy  spread  from  the  source  which  operates  at  about  12  00°C.     Mass 
analysis  shows  that  beam  purities  are  well  in  e 
isotopically  enriched  emitter  material  is  used. 


7     + 
analysis  shows  that  beam  purities  are  well  in  excess  of  99%     Li    when 


II.     BEAM  DYNAMICS 

The  beam  geometry  in  the  absence  of  space  charge  effects  is  shown 

in  Fig.  2.    As  the  beam  leaves  the  aperture  of  radius  R    the  rays  ema- 

B 

nating  from  the  edge  are  initially  directed  as  shown  and  those  from 

interior  points  are  contained  within  the  outer  rays  so  as  to  be  focused 

toward  a  spot  of  radius  r    a  distance  L  away.    The  origin  of  coordinates 

is  such  that  the  beam  travels  along  the  z  axis,  the  coordinated  of  the 

aperture  is  z     and  of  the  image  z     .    The  beam  is  assumed  to  be  axially 
B  a 

symmetric. 

In  considering  the  space  charge  effects  we  shall  follow  the  work 
of  Glaser    and  the  analyses  of  Watson    and  Wendt     .    We  shall  assume 
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FIGURE    2 


a  beam  of  charged  particles  of  sufficient  density  to  be  considered  a 
continuous  medium  for  purposes  of  constructing  energy  functions. 

The  effect  of  the  space  charge  of  the  beam  is  to  deflect  the  rays 
away  from  the  axis.    We  assume  that  the  current  density  is  uniform 
across  the  beam  cross  section  and  zero  outside  the  outer  rays,  and  that 
the  total  current  is  constant.     If  the  convergence  of  the  rays  is  suffi- 
ciently small  then  a  small  segment  of  the  beam  forms  a  cylinder  of 
nearly  constant  radius  in  which  the  axial  forces  are  nearly  zero,  and 
hence  the  axial  velocity  z  is  very  nearly  constant  and  given  by 
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T-eSV 


(1) 


where  e  =  charge  on  the  particle,  V  =  beam  potential,  and  m  =  particle 
mass . 

We  shall  now  proceed  to  calculate  the  trajectories  of  the  outer 
and  inner  rays  as  shown  in  Fig.  2  under  the  influence  of  space  charge 
(in  field  free  space)  subject  to  the  simplifying  assumptions  given  above. 

Trajectory  of  the  Outer  Ray 

The  radial  component  of  the  electric  field  due  to  the  space  charge 
can  be  written  .. 

where  i  is  the  total  current,  eV  is  the  beam  energy,  and  r     is  the  radius 
of  the  outer  ray  at  any  point.     Let 

*\_  *    "vvT 

be  the  charge  to  mass  ratio  of  the  particles.    The  radial  equation  of 

motion  of  a  particle  on  the  outer  ray  is  then 


Wr>     *h*FT7       77==  (3) 


If  we  define  the  constant  K.  to  be 


t  K.  to  ^ 

1       « 

K,  =   -4, 


we  may  write  Eq.    (3)  as 
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r*'*^ 


Multiplying  both  sides  by  dr    and  integrating  gives 
where  the  constant  of  integration  r    is  chosen  such  that  at 

O 

Using  equations  (1)  and  (4)  we  can  write         i 

where  we  choose  the  negative  root  in  order  to  restrict  consideration  to  the 
portion  of  the  beam  between  the  aperture  and  the  focal  point,  i.e.  ,  the 
region  where  the  beam  is  convergent  and  r    <  0.    This  gives  on  integration 

h      r-  .      _     V 
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we  find 


6  Ju  Lv  **-  WJ    ^ 

where  the  primes  denote  variables  of  integration. 
Finally,  setting  I 

,4 

*-*_  =  -  2K.        &*p  \  s    \  oU; 

where  8  Jf  ^      L  J 

If  we  define  the  Dawson  function        ' 

DCs)-    Jo   ^{s'1}^' 

then  Eq.   (9)  becomes 


(8) 
(9) 
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(10) 


(11) 


Since  tables  of  D(S)  are  available       we  can  evaluate  equation  11  in 

terms  of  the  constant  K.     This  equation  gives  r    as  an  implicit  function 

of  z  and  the  parameters  R_,  r  ,L,  wand  the  beam  perveance  P. 

B      a  v 

To  express  K  in  terms  of  the  minimum  radius  r  o  and  the  beam 
perveance  "P  9  rr       we  note  from  above,  that 

K  -  v;  {^~  (12) 


where  K  ,    - 


Thus 


or 


U  (4) 

,>»/  , (13) 

K-O.P      *  where     Q     „   yfi^J^ 

Trajectory  of  the  Inner  Ray 

We  know  from  Gauss'  Law  that  only  the  charge  enclosed  by  the 
inner  ray  will  affect  its  trajectory.    Assuming  a  uniform  radial  charge 

distribution,  the  ratio  of  the  current  included  within  r9  to  the  current 

2 
included  within  r,  will  be    r2         and  since  the  field  depends  only  on  the 

enclosed  charge  we  can  use  Eqs.   (2)  and  (4)  to  write 

and  the  radial  equation  of  motion  of  a  particle  on  the  inner  ray  is 

£  -  1**7*  (15) 

where  r^  is  given  by 

from  Eq .   (8) . 

We  may  differentiate  r«  parametrically  as 

dub  "  "cTF  dr,  "aa^  ■  "st"*  **t  ' 

where  £    is  given  by  (6).    ds/dr    may  be  found  from  (8),  giving 


r  -  (  i — ")  e**  L£i  *k 


and 
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(17) 


Putting  (18)  in  (15)  gives 

f£%-  2s^X±  -2ru  -O  (19) 

which  has  the  solution^* 
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where  C-^  and  C2  are  constants  of  integration,  and 

S 


0  (s)  is  the  well-known  error  function  which  is  extensively  tabulated   ^  . 

We  evaluate  the  constants  of  integration  by  noting  that  at  z  =  z 
(see  Fig.  2) 

and 
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Note  also  that  from  equation  (9)  j         t^ 

Applying  Leibnitz' rule  gives  «^0  *-        J 
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Applying  Leibnitz1  rule  to  (21)  gives 
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and  . 

Putting  (23),   (28),   (29),  and  the  second  equation  of  (22)  into  (27)  and 
evaluating  at  ^  =  Rg  ,  gives 


5£i.  £[«.(. ♦C.O^e,.,  {- 


sl 


where  0B  denotes  0  (Sg ) .  Solving  (2  6)  for  C-i  and  putting  this  and  (24)  in 
(30)  gives  . 

1  +  Cj.^g    "  To1   L 

™  v.*L  (31) 

With  this  we  have  immediately 

C2    =      ^-j—j (32) 

and  from  (26)  rS 

so  the  complete  solution  for  r,  (see  equation  20)  is  given  by 


^-n.txpi^lj-N^-^)] 


(33) 
(34) 


where  r    is  given  by  (25),  N  by  (31),  s_  by  (24)  and   <t>  by  (21).     It  should 
o  B 

be  noted  that  equation  (34)  gives  rasa  function  of  the  constants  R    , 

&  B 

r    ,  L,  tf  and  the  beam  perveance  P  at  every  value  of  z. 
a 


Calculation  of  Minimum  Beam  Radius  and  Postion 


The  spreading  of  the  rays  under  the  influence  of  space  charge  is 

illustrated  in  Fig.  3.    We  shall  distinguish  between  two  cases  depending 

on  whether  the  value  of  the  beam  perveance  P  is  larger  than  or  smaller 

than  a  critical  perveance  P      defined  below: 

cr 

Case  1;    P>P      (Fig.  3a):  For  this  case  the  outer  ray  reaches 

its  minimum  before  its  intersection  with  the  inner  ray 

originating  at  the  diametrically  opposite  side  of  the 

aperture.     Thus  the  outer  ray  reaches  its  minimum 

value  r    =  r    before  -r„  =  r.  and  hence  the  minimum 
1        o  2        1 

radius  occurs  at  z  =  z    . 

o 

Let  z    -z„  =  L„  .    At  rn  =r    ,  s  =  0  (from  equation  8). 

0       B         F  1        ° 

Since  D(o)  =  0  we  get  from  equation  (11) 

Lf     '     i*<  £>(sb}  (35) 

and  ro  =  r     is  given  by  equation  (25). 

Case  2;    P<P      (Fig.  3b):    Here  the  outer  ray  will  be  intersected 
cr  

by  the  inner  ray  originating  at  the  diametrically  opposite 
side  of  the  aperture,  before  r    reaches  its  minimum  value 
r  .     If  we  call  this  point  of  intersection  z     then  at 

Z  =  ZF'--r2=ri=V 

Let  z     -  z     =  L     .    The  minimum  radius  occurs  where 

r]  =   -v   .    Equation  (16)  and  (34)  gives 

from  which 

2__  (36) 
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FIGURE    3 
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s  can  be  found  from  (36)  and  we  find  r„  from 

F  F 

rp     =    rQ   exp    £  sp2|  (37) 

and 

Lp     =    2  K    [D  (sb)  -  p  (sp)]  (38) 

When  P  =  P       ,  Lp  (case  1)  =  Lr,  (case  2)  and  the  method 
cr         F  r 

of  either  case  may  be  employed  in  evaluating  rp  and  L     . 

III.    RESULTS  OF  NUMERICAL  COMPUTATIONS 
In  this  section  we  present  the  results  of  numerical  computations  for 
the  image  distance  L     (Equ.  35  and  38)  and  the  image  radius  r„  (Equ.  25 
and  37).     Each  of  the  figures  (Fig.  4-18)  shows  3  curves  of  Lp  as  a 
function  of  beam  perveance       r        ( ~     "1_  Up    J  for  three  values  of 

the  beam  convergence  angle  y  (Fig .  2)  (tan  y-  0.  03,  0  .09  and  0.15). 
On  the  same  figure  we  plot  r     as  a  function  of  P  for  the  same  three 
values  of  tan  y .    The  curves  are  aligned  so  that  values  of  Lp  and  rp 
can  be  read  for  any  given  value  of  P  and  the  same  value  ©f  tan  y  by 
drawing  a  vertical  line  at  the  proper  value  of  perveance.    Three  values 
of  the  final  aperture  radius  R„  (See  Fig.  2)  (Rg  =  0.1,  0.25  and  0.50 
cm)  and  five  mass  values  were  used:    electrons,     H    ,      Li    ,  "   Ne 
and       Rb    . 
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